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We study ground-state properties of the two-dimensional random-bond Ising model with couplings 
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700. We study the behavior of the domain-wall energies and of the magnetization. We find that the 
paramagnet-ferromagnet transition occurs at pc '~ 0.103 compared to the concentration p„ 0.109 
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I. INTRODUCTION 

Despite more than two decades of intensive research, 
many properties of spin glasses^, especially in finite 
dimensions, are still not well understood. For two- 
dimensional Ising spin glasses it is now widely ac- 
cepted that no ordered phase for finite temperatures 
gjj.jg^g2j3j4j5j6_ Furthermore, it seems clear that the behav- 
ior can be described well by a zero-temperature droplet 
scaling approac h but one needs quite large sys- 
tem sizes to observeii the true behavior. One unan- 
swered question is whether an additional phase, usually 
called random antiphase, exist3i2iiiilii^ for T in two 
dimensions with an asymmetric distribution of random 
bonds. Also it was no clear whether at low tempera- 
tures the phase boundary of the ferromagnetic phase is 
perpendicular^^' to the p axis, p denoting the concen- 
tration of the antiferromagnetic bonds. The aim of the 
present paper is to reinvestigate this issues by studying 
the domain-wall energy and magnetization at zero tem- 
perature via the determination of exact ground statesi^ 
for large system sizes and huge sample numbers. This 
allows us to draw much more reliable conclusions in com- 
parison to past studies, where only considerable smaller 
system sizes could be studied. 

The organization of the paper is as follows. In the next 
section, we will expose the model and briefly describe the 
polynomial matching algorithm, which allows us to treat 
large system sizes. Then we will present our results for 
the domain-wall energy. In the succeeding section, we 
explain the additional methods used to obtain the mag- 
netization and show the results. Finally, we summarize 
and draw our conclusions. 



II. THE MODEL AND THE METHOD 

The model consists oi N — spins Si = ±1 on a 
simple square lattice with periodic boundary conditions 



direction. The Hamiltonian is 

H = — "^^JijSiSj (1) 

where the sum runs over all pairs of nearest neighbors {ij) 
and the Jy are the quenched random spin-spin couplings. 
The couplings are set independently antiferromagnetic 
{ Jij — —1) with a probability p e [0, 1] and ferromagnetic 
{Jij = +1) with probability (1 —p). 
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FIG. 1: The phase diagram of the two-dimensional random- 
bond Ising model, with the concentration p of antiferromag- 
netic bond and the temperature T on the vertical axis. It has 
been conjectured that the phase boundary from the NP falls 
vertically to the p-axis. 

The phase diagram of the model as a function of tem- 
perature T and the concentration p is shown in Fig. ^ 
The pure system {p — 0) has a transition at a Curie tem- 
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paramagnetic. When antiferromagnetic bonds are intro- 
duced {p > 0), the ferromagnetic phase is destroyed at 
a threshold concentration p^iT). A particular curve on 
the p — T plane is known as the Nishimori linei^ (^^L). 
It is defined by the equation exp2/3 = On this line 
the internal energy is analytic and the spin-spin correla- 
tion functions obey the equalities {aiGj)"^^ ^ = (aiaj)'^'^ , 
for integer k. It was also provenSfl that a multicritical 
point delimiting two critical behaviors on the ferro-para 
boundary coincides with the intersection of the NL with 
the boundary. Besides, by studying domain- wall energies 
of exact ground states for system's sizes up to L = 32, 
Kawashima and Rieger^^ found that the stability of the 
ferromagnetic and spin-glass order cease to exist at a 
unique concentration for the antiferromagnetic bonds, so 
they concluded that there is no intermediate spin-glass 
phase. 

In this paper we want to compute numerically with 
high accuracy the critical concentration pi^^ = pc{T ^ 0) 
corresponding to the para-ferro transition at zero temper- 
ature, and to compare this result with the believed value 
of the Nishimori point p„ = 0.109(1)^^1^^. Furthermore, 
we want to check with high accuracy, whether there is an 
intermediate spin-glass phase at nonzero temperature. 

We can reach a much higher precision compared to pre- 
vious studies, by applying a matching algorithm. This 
allows to compute exact ground states for large system 
sizes, N = 700^ spins in our case. Let us now explain 
just the basic idea of the matching algorithm, for the 
details, see Refs. .23,,„24,25. The method works for spin 
glasses which are planar graphs, this is the reason, why 
we apply periodic boundary conditions only in one direc- 
tion. In the left part of Fig. [21 a small 2d system with 
open boundary conditions is shown. All spins are as- 
sumed to be "up" , hence all antiferromagnetic bonds are 
not satisfied. If one draws a dotted line perpendicular 
to all unsatisfied bonds, one ends up with the situation 
shown in the figure: all dotted lines start or end at frus- 
trated plaquettes and each frustrated plaquette is con- 
nected to exactly one other frustrated plaquette. Each 
pair of plaquettes is then said to be matched. Now, one 
can consider the frustrated plaquettes as the vertices and 
all possible pairs of connections as the edges of a (dual) 
graph. The dotted lines are selected from the edges con- 
necting the vertices and called a perfect matching, since 
all plaquettes are matched. One can assign weights to 
the edges in the dual graph, the weights are equal to 
the sum of the absolute values of the bonds crossed by 
the dotted lines. The weight A of the matching is de- 
fined as the sum of the weights of the edges contained in 
the matching. As we have seen, A measures the broken 
bonds, hence, the energy of the configuration is given by 
E = j) I'^ijl +2A. Note that this holds for any con- 

figuration of the spins, since a corresponding matching 
always exists. Obtaining a ground state means minimiz- 
ing the total weight of the broken bonds (see right panel 



matching. This problem is solvable in polynomial time. 
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FIG. 2: 2d spin glass with all spins up (left, up spins not 
shown). Straight lines are ferromagnetic, jagged lines antifer- 
romagnetic bonds. The dotted lines connect frustrated pla- 
quettes (crosses). The bonds crossed by the dotted lines are 
unsatisfied. In the right part the GS with three spins point- 
ing down (all other up) is shown, corresponding to a minimum 
number of unsatisfied bonds. 

The algorithms for minimum-weight perfect 
matchingsSSiSi are among the most complicated al- 
gorithms for polynomial problems. Fortunately the 
LEDA library offers a very efficient implcmcntation^iSi, 
except that it consumes a lot of memory, which limits in 
our case the size of the systems to about N — 700^ on a 
typical 500 MB workstation. 



III. DOMAIN- WALL ENERGIES AND 
FINITE-SIZE SCALING. 



We calculate the domain-wall energy 6E defined by 
SE = Ep~ Ea where Ep and Ea are the ground-state en- 
ergies with periodic and the anti-periodic boundary con- 
ditions in the a:-direction, respectively. We take an aver- 
age over the disorder. We are interested in the exponents 
p and 9s that characterize the system-size dependence of 
the mean AE and the width a{6E) of the distribution of 
the domain- wall energies: 



AE cx LP and a{SE) cx 



(2) 



For a general dimensions d of the system, a positive value 
of p indicates the stability of a ferromagnetic phase For 
p < 0, no ferromagnetic ordering is present. Then, in 
dimension d above the lower critical dimension dc, we 
have 9s > and spin glass-ordering is stable against 
thermal fluctuations. On the contrary, when ^5 < 0, 
thermal fluctuations prevent spin-glass orderinglfl. The 
current believe is that in d — 2, 9s < holds for all 
concentrations pc^'' < p < 1 — pi^"^ . 

We have computed a{5E) and AE for sizes up to 
L = 700 and for values of p ranging from 0.100 to 0.107. 
We performed an disorder average of a number of re- 
alizations ranging from 30000 for the smallest sizes to 
typically 2000 for the largest size L — 700. In Fig. Oand 
0]the mean and the width of the distribution of domain- 
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FIG. 3: The domain-wall energy AE of the bond-random 
model plotted as a function of the system size L for various 
ferromagnetic-bond concentrations p. 
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FIG. 4: The width of the distribution of the domain-wall en- 
ergy a{5E) of the bond-random model plotted as a function 
of the system size L for various ferromagnetic-bond concen- 
tration p. 



We denote by pi^^ and pc^'' the critical concentrations of 
antiferromagnetic bonds at which the asymptotic L de- 
pendences of AE and (j(6E), respectively, change from 
increasing to decreasing, i.e., the concentrations where a 
ferromagnetic phase and a spin-glass phase, respectively, 
cease to exist at finite temperature. We conclude from 
the figures that pi^^ ~ 0.103, while for pi?'' the "tran- 
sition" is less sharp but the value is between 0.103 and 
0.105. For small sizes, the width even seems to increase 
for all values of p € [0.1,0.107] we have considered. For 
small sizes, at intermediate concentrations p € [0.1,0.15], 
the mean of the domain-wall energy already decreases 
with system size, while the width of the distribution first 
increases, which appears as if the system exhibits some 
kind of spin-glass phase. This is probably the reason that 
in some previous studiesiSiiii4 the existence of an addi- 
tional intermediate phase has been assumed. We see that 
we have to consider large system sizes to observe the true 
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FIG. 5: The scaling plot of mean domain-wall energy AE 
versus the concentration of the antiferromagnetic bonds, using 



the values p, 
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FIG. 6: The scaling plot of a{SE) versus the concentration 
of the antiferromagnetic bonds, using the values p, 
(j>2 = 0.74, ^2 = -0.13. 



(2) 



0.104, 



Another way to compute pi^^^^^'' is to check the scaling 
relations for AE and (t{SE) proposed ini^: 



a{SE)L'^^=M{p~-p(?y}L'^^] 



(3) 



(4) 



The parameters Pc, 4> and ■0 for both moments of do- 
main wall energies have to be chosen such that a good 
data collapse for all data is obtained. To quantify the 
"goodness" of this fit, we used an appropriate cost func- 
tion S{pc,(j),4') introduced in whose minimum value 
should be close to unity when the fit is statistically ac- 
ceptable. To minimize S'(pcj0, "0) we used the imple- 
mentation of the simplex method offered by Numerical 
Recipes librarji22,. The best fits give the estimates 
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with S 0.75 and 

^(2) ^ 0.104(2), (j)2 = 0.74(5), ^2 = -0.13(5) (6) 

with S = 0.65. The resulting scahng plots are shown in 
Figs. El and O We have estimated the error bars given 
above in the following way. For each parameter, we fix it 
to different values and perform the minimization over the 
remaining two parameters. In Fig. [7|we show as example 
a plot of this partly minimized value of S{pc, ip, (f>) as a 
function of pc- Our error bars are the ranges of values 
where S{pc, ip, (j)) increases to twice of its minimum value. 
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FIG. 7: Plot of the minimum value of S{pc, tp, 0) for different 
fixed values of pc- 

Within the statistical errors the critical parameters for 
both moment of SE agree: this strongly suggest the ab- 
sence of a spin-glass phase. Therefore there is a discrep- 
ancy between the critical concentration p evaluated at 
the NP and at zero temperature, in disagreement with 
the conjecture p„ — p^c^^ by Nishimori''-": and later by 
Kitatan:^^. Later, Le Doussal and Harris LeDoussalHar- 
ris have shown that the tangent to the phase boundary 
at the NP is vertical. But this result does not exclude 
the possibility of a reentrance in the phase diagram as 
shown in Fig. ^ 



The problem one has to face when studying system 
with discrete interactions, is the exponentially large num- 
ber of states all giving the same energy. Hence there is 
no unique ground-state magnetization. For a given set 
of bonds we here determine one exact ground state using 
an efficient polynomial time "matching" algorithm, but 
we are not able to enumerate all the ground states^. 

In order to check if "typical" configurations with re- 
spect to the magnetization are found, we first performed 
a zero-temperature Monte Carlo (MC) simulation which 
consists in flipping all spins with zero local field, starting 
with the ground-state configuration. This allows to ex- 
plore all states within a single-spin-flip cluster of ground- 
state configurations (see below) . In Fig |S1 a typical evo- 
lution of the magnetization, here for L — 300, p = 0.1, 
averaged over 100 samples, is shown as a function of the 
number of Monte Carlo steps. We observe that after few 
hundred MC steps the simulation is equilibrated. Addi- 
tionally, we see that the value of the magnetization found 
by the ground-state algorithm, i.e. for zero MC steps, is 
slightly outside the scattering of the datapoints for large 
Monte Carlo steps. But compared to the final statisti- 
cal error bars (see below), this difference is negligible. 
Hence we conclude that the ground state obtained by 
the matching algorithm exhibits a typical magnetization 
of the cluster, in which the ground state is located. 
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IV. STUDY OF THE MAGNETIZATION. 

We furthermore study the para-ferro transition by eval- 
uating the magnetization and using the Binder cumulant 
crossing method^ii^^. The Binder cumulant is given by 

B(p,L) = - 1--^ , (7) 



where m — j^J^'^i is the magnetization and (. . .) de- 
notes the average over the disorder. For second order 
phase transitions, the curves for different sizes intersect 
at one point, the critical concentration pi^\ This is a 

_ri2^-ij__ 1- 



FIG. 8: Magnetization averaged over 100 samples {L — 300, 
p = 0.1) as a function of the number of MC steps for a zero- 
temperature single-spin-flip dynamics. 

Anyway, the set of ground states usually is divided into 
several clusters^i: different ground states belong to the 
same cluster if they are related by a sequence of single 
spin flips, each leaves the energy the same. Ground states 
in different cluster can only be reached from each other by 
making cooperative flips of multiple spins or when using 
single-spin flips via increasing the energy. This means 
that with our single-spin-flip MC algorithm at T = 
the system always stays in the same cluster. In prin- 
ciple one can enumerate all ground statesSi. Since the 
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the system size, only small systems can be treated like 
this. 

Thus, we have applied an alternative method to find 
different ground states, the so called e-coupling method. 
It allows to obtain ground states from different clusters, 
but no exhaustive enumeration is necessary. The basic 
idea is to first add a perturbation to the system which 
tends to increase the energy if two neighbouring spins 
are in the same relative orientation like in the ground 

state and then to recalculate the ground state. Let S'j^'' 
be the ground-state spin configuration. We then perturb 
the couplings Jij by an amount proportional to S^^^ S^^^ 
in order to repel the system from the ground state. This 
perturbation, which depends upon a positive parameter 
e, is defined by 

^ J:,^- + A J„ (8) 

where 

AJ.,=-^S^^S^ (9) 

where Nb is the number of bonds in the system. We 
then recompute the ground state and check that the new 
configuration is still a ground state of the unperturbed 

(2) 

Hamiltonian. This is our second ground state . The 
next step, to obtain a third ground state, consists in 
adding a perturbation in order to repel the system from 
both ground states obtained so far. This process can be 
iterated. For a number n of steps of this process, we 
have: 

n 

^J^^-'^Ew/^'^sf, (10) 

k—l 

where a is a scaling factor that we choose equal to 1. 
In this way we hope to find configurations belonging 
to different clusters, even if this procedure is not com- 
pletely under control since it is a biased random sam- 
pling in the space of configurations. To test the behav- 
ior of our method, we have calculated for each step k 
of the e-coupling approach the magnetization m'^ of the 
A;th ground state. Fig. |51 shows the resulting behavior of 
[m*^], [. . .] denoting the average over the first k—l itera- 
tions of the e-coupling method and over 1000 configura- 
tions the disorder {L — 300, p = 0.1). We observe that, 
within the statistical error bars, the magnetization for 
the first ground state, indicated by the horizontal lines, 
agrees well with the value obtained after averaging over 
several different degenerate ground-state configurations. 
Note that for small k the difference is larger. This is due 
to the fact that the e-coupling method repels each con- 
figuration from the previously obtained ground states. 
Hence, for small k it will move strongly away from the 
first ground state. With increasing number k of steps, 
the different ground states will be scattered around in 
configuration space. If the k = 1 ground state is typi- 



will converge to the initial value again, which seems to 
be the case here. Note that we have also checked for the 
different iterations of the e-coupling method that per- 
forming an additional T = MC simulation changes the 
obtained magnetization values again only slightly. Since 
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FIG. 9; Magnetization obtained after averaging over k inde- 
pendent ground states as a function of k. The horizontal lines 
indicate the magnetization of the first ground state k — l. 
The data is for L — 100, p — 0.103 and averaged over 1000 
realizations. 

the e-coupling method strongly slows the simulations and 
the differences are negligible within statistical error bars, 
we have restricted the simulations to the immediately ob- 
tained ground states (fc = 1), i.e. without applying the 
e-coupling method. 

In Fig. 1101 the Binder cumulant is shown as a function 
of the concentration of p for different sizes L < 500, where 
we have obtained data for all concentrations p — 0.1, 
0.101, . . . , 0.107. It intersects close ~ 0.103 expect for 
the largest size, where the stistics is not so good. Hence, 
we again conclude p]. = 0.103(1) which agrees well with 
what we have obtained above by studying domain wall 
energies. This also indicates that the matching algorithm 
indeed finds ground states which are typical with respect 
to the magnetization. 

Finally, we present further ways to check the previous 
conclusions. We have performed another treatment of 
the data by trying to collapse all curves in Fig. ^| in a 
single one. This means we want to find the parameters 
which satisfy the finite-size scaling relation for the Binder 
cumulant: 

B{p,L)^B{L^'^{p-p,)) (11) 

where v is the critical exponent of the correlation length. 
We vary pc and v in order to minimize the functional 
5(pc, i^)- We find its minimum value for pc — 0.103 and 
V — 1.55(1) {S ^ 2.2), the resulting data collapse is illus- 
trated in Fig. ^] The value for v is consistent with the 
value V = 1.50(3) at the NP reported in Ref. HI but dif- 
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FIG. 10: Binder cumulant B{p,L) = |(1 - {m)"^/{m'^f) as a 
function of the concentration of antiferromagnetic bonds p. 



in the work of Merz and Chalkcr much larger system sizes 
are treated compared to the work of Honecker et. al., the 
value v = 1.50(3) appears more reliable. Hence, our re- 
sult indicates that the transition at NP and at T = 
are in the same universality class. 
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FIG. 11; Scaling plot of the Binder cumulant B(p,L) as a 
function of (p — pc)L^^" with pc — 0.103 and u = 1.55. 
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FIG. 12; Scaling plot of the rescaled magnetization mL^^'^ as 
a function of {p-pc)L'^^'' with pc = 0.103, v = 1.55, /3 = 0.9. 



By using a matching algorithm, we could study systems 
which are much larger than in previous studies. 

We find that both ferromagnetic and spin-glass phases 
cease to exist at the same concentration pc — 0.103(1). 
This means, we do not find any sign for an intermedi- 
ate "random antiphase" . Furthermore, the values of Pc 
at T = and at the Nishimory point, found in the most 
reliable studies so fariSLSS, are different, indicating a reen- 
trance of the paramagnetic phase. 

Large system sizes L > 500 are needed to observe 
the true thermodynamic behavior with good accuracy. 
Slightly above pc, the width of the distribution of 
domain- wall energies increases for small sizes, while it 
decreases for larger sizes. Note that in principle we can- 
not exclude that a similar turnover happens for smaller 
concentrations, e.g. p = 0.103, at even larger system 
sizes L > 700, which are out of reach of our algorithm. 
Nevertheless, this would mean that the real Pc is even 
smaller, hence the discrepancy between pc and p„ would 
increase and the reentrance became stronger. 



Finally, we study the finite-size scaling behavior of the 
magnetization. The prediction for the magnetization is 

m{p) = L^f^^^rndp - Pc)L^^'') (12) 

where /3 is the critical exponent of the magnetization. 
Using Pc = 0.1032, i/ = 1.55, a good data collapse is 
obtained with /3 = 0.9(1), see Fig. [HI 

V. CONCLUSIONS 

To summarize, we have performed a systematic calcu- 
lation of the domain-wall energy at zero temperature for 
the two-dimensional random-bond Ising model for dif- 
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